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' Symmetry and its spontaneous breaking is a central theme in modern physics. Perhaps no symmetry is more 

■ fundamental than time translation symmetry, since time translation symmetry underlies both the reproducibility of 
D ] experience and, within the standard dynamical frameworks, the conservation of energy. So it is natural to consider 

■ the question, whether time translation symmetry might be spontaneously broken in a closed quantum-mechanical 

■ system. That is the question we will consider, and answer affirmatively, here. 

y—{ . Spontaneous breaking of spatial translation symmetry is of course very common, typically leading to crystals. Here 
we are considering the possibility of time crystals. Classical (and semiclassical) time crystals are considered in a 
companion paper in this note the primary emphasis is on quantum theory. 

Several considerations might seem to make the possibility of quantum time crystals implausible. The Heisenberg 
equation of motion for an operator with no intrinsic time dependence reads 

where the last step applies to any eigenstate oi H. This seems to preclude the possibility of an order parameter 
that could indicate the spontaneous breaking of infinitesimal time translation symmetry. Also, the very concept of 
"ground state" implies state of lowest energy; but in any state of definite energy (it seems) the Hamiltonian must act 
^ ■ trivially. Finally, a system with spontaneous breaking of time translation symmetry in its ground state must have 
^\ \ some sort of motion in its ground state, and is therefore perilously close to fitting the definition of a perpetual motion 
CO ■ machine. 

in 

(N 

^ ■ I. Ring Particle Model 

o: 

■ And yet there is a familiar physical phenomenon that almost does the job. A superconductor, in the right circum- 
stances, can support a stable current-carrying ground state. Specifically, this occurs if we have a superconducting 
ring threaded by a flux that is a fraction of the fiux quantum. If the current is constant then nothing changes in time, 
so time translation symmetry is not broken; but clearly there is a sense in which something is moving. 

We can display the essence of this situation in a simple model, that displays its formal structure clearly. Consider a 
] particle with charge q and unit mass, confined to a ring of unit radius that is threaded by fiux 2'Ka/q. The Lagrangian, 
canonical (angular) momentum, and Hamiltonian for this system are respectively 

1 •„ 

L = —(f) + aip 
TT^ = 4) + a 

H = \{^^~af (2) 

7r0, through its role as generator of (angular) translations, and in view of the Heisenberg commutation relations, is 
realized as —i^- Its eigenvalues are integers I, associated with the states \l) — e*'"^. For these states we have 

= I- a 

{mi) = (3) 
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The lowest energy state will occur for the integer Iq that makes I — a smallest. If a is not an integer, we will have 

(^0 101^0) = lo-a ^ (4) 

The case when a is half an odd integer requires special consideration. In that case we will have two distinct 
states |a ± i) with the minimum energy. We can clarify the meaning of that degeneracy by combining two simple 
observations. First, that the combined operation Gk of multiplying wave functions by e*'^'^ and changing a a + k, 
for integer fc, in the Lagrangian leaves the dynamics invariant. Indeed, if we interpret a in i as embodying a constant 
gauge potential, Gk is a topologically non-trivial gauge transformation on the ring, corresponding to the multiply- 
valued gauge function A ^ A + VA, A = k(j)/q. Note that the total flux is not invariant under this topologically 
non-trivial gauge transformation, which cannot be extended smoothly off the ring, so L is modified. Second, that the 
operation of time- reversal T, implemented by complex conjugation of wave- functions, takes \l) — | — /) and leaves 
the dynamics invariant if simultaneously a — )■ —a. Putting these observations together, we see that the combined 
operation 

f - GaaT (5) 

leaves the Lagrangian invariant; it is a symmetry of the dynamics and maps |^) — >■ \2a — l).T interchanges |^ ± a) — > 
1^ =F a). Thus the degeneracy between those states is a consequence of a modified time-reversal symmetry. We can 
choose combinations |a-fi)±|a — -i) that simultaneously diagonalize H and T; for these combinations the expectation 
value of 4> vanishes. 

Returning to the generic case: For a that are not half-integral time-reversal symmetry is not merely modified, but 
simply broken, and there is no degeneracy. How do we reconcile (ZqI^I^o) 7^ with Eqn. ([l}? The point is that 0, 
despite appearances, is neither the time-derivative of a legitimate operator nor the commutator of the Hamiltonian 
with one, since 0, acting on wave functions in Hilbert space, is multivalued. By way of contrast operators corresponding 
to single- valued functions of 0, spanned by trigonometric functions Ok = e**^*, do satisfy Eqn. ^ for the eigenstates 

I*) = 10- 

Wave functions of the quantized ring particle model correspond to the (classical) wave functions that appear in the 
Landau-Ginzburg theory of superconductivity. Those wave functions, in turn, heuristically describe the wave function 
for macroscopic occupation of the single-particle quantum state appropriate to a Cooper pair, regarded as a particle. 
Under this correspondence, the non-vanishing expectation value of (j) for the ground state of the ring particle subject 
to fractional flux maps onto the persistent current in a superconducting ring. 



II. Symmetry Breaking and Observability 

As mentioned previously, the choice of a ground state that violates time translation symmetry r must be based 
on some criterion other than energy minimization. But what might seem to be a special difficulty with breaking r, 
because of its connection to the Hamiltonian, actually arises in only slightly different form for all cases of spontaneous 
symmetry breaking. Consider for example the breaking of number (or dually, phase) symmetry. We characterize such 
breaking through a complex order parameter, that acquires a non-zero expectation value, which we can take to be 
real: 

(0|i>|0) = w 7^0 (6) 

We also have states jcr) related to |0) by the symmetry operation. These are all energetically degenerate and 
mutually orthogonal (see below), and satisfy 

(a|$|a) = ve'" (7) 

The superposition 

l^^> - ^ / dcr \<t) (8) 



is energetically degenerate with all the |(t), and it is symmetric, with 

{n\<^\n) = (9) 
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Why then do we prefer one of the states \a) as a description of the physical situation? The reason for that preference 
is closely related to reason for the orthogonality of the different \a) states, which we now recall. We envisage that 
our system extends over a large number N of identical subsystems having correlated values of the long-range order 
parameter cr, but otherwise essentially uncorrelated. Then we can express the total wave function in the form 

N 

^a{xi,...,XN) ~ J|V<T(a;j) (10) 

For different values ct, a' we have therefore 

N 

(*.'|*.) « X{{^.'{xi)\i^.{xj)) = ifa'.f->0 {a'^a) (11) 

for large N, since Ifa'al < 1- Similarly, for any finite set of local observables (that is, observables whose arguments 
include only upon a finite subset of the Xj), we have 

Since the off-diagonal matrix elements vanish, any world of local observations (including "observations" by the envi- 
ronment) can be described using a single \a) state. Averaging over them, to produce is a purely formal operation. 
Measurement of a non-singlet observable will project onto a \a) state. 

This analysis [2] brings out several relevant points. The physical criterion that identifies useful "ground states" is 
not simply energy, but also robust observability. Mathematically, that requirement is reflected in the orthogonality 
of the Hilbert spaces built upon \a) states by the action of physical observables. The large N limit is crucial for 
spontaneous symmetry breaking. It is only in that infinite degree-of-freedom, or as it is usually called infinite volume, 
limit, that the |j) and |cr) states are all degenerate, and the latter are preferred. And of course the overarching 
point is that the preceding discussion applies, with only symbolic changes, when we consider possible breaking of 
time-translation r in place of phase symmetry. 



III. Soliton Model 

After these preparations, it is not difficult to construct an appropriate model. We consider a large number of 
ring-particles with an attractive interaction. Heuristically, we can expect that they will want to form a lump and, in 
view of Eqn. (j?]), that they will want to move. So we can expect that the physical ground state features a moving 
lump, which manifestly breaks r. 

In order to make contact with the argument of the previous section, we should have an appropriate notion of locality. 
For simplicity let us assume that the particles have an additional integer label, besides the common angle 0, and that 
the physical observables are of finite range in the additional label. (We might imagine an array of separate rings, 
displaced along an axis, so that the coordinates of particle j are {(f>,x — ja).) We do not insist that the Hamiltonian 
is local in that sense, however. I will return to these conceptual issues below, after describing the construction. 

An appropriate Hamiltonian is 

N X ^ ^1 

with the understanding that H acts on periodic functions, so the 6 interaction is well defined. (Here the discrete 
index appears as a subscript.) 

We work in the mean field approximation, taking a product ansatz 

N 

vI/(^i,...,</,„) = Y[ij{cPj) (14) 
and solving an approximate one-body equation for i/j. To get such an equation, we define an effective potential 

N 

W{cl)j) = fl[d<l>kr{MVi^{M (15) 
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so that 



(^'IKff.l'J') - (*|^^|*) (16) 

Then the effective Schrodinger equation for 5*, 

f)\[ ^ 

.7 = 1 



reduces to the one-body non-hnear Schrodinger equation 

'dt ~ 2 



= i(^^-a)2^-A|V|V (18) 



for ip. 

Consider first the case a = 0. Eqn. (IT51) can be solved for a stationary state in terms of the Jacobi dn ehiptic 
function, with 

V'o(0) = r dn (r\/A0, fc^) 

f = -r2A(l-:^) (19) 

with P a disposable parameter. To fix the parameters fc, r we must impose 27r periodicity in (j) and normalize "00 ■ 
Those conditions become 



K{k^) ^ nr^/X (20) 



in terms of the complete elliptic integrals E{k^),K{k^). We can solve E{k^)K{k^) = ^ for k^, given A. The 
minimum value of the left-hand side occurs at fc = and corresponds to A = ^. Here dn {u, 0) reduces to a constant, 
and £ — —1/4. As A increases beyond that value k rapidly approaches 1, as does E{k^). dn(u, fc^) — )■ sechu and 
£ — >■ — A^/8 in that limit. Of course the constant solution with £ — —X/2tt exists for any value of A, but when 
A exceeds the critical value the inhomogeneous solution is more favorable energetically. These results have simple 
qualitative interpretations. The hyperbolic secant is the famous soliton of the non-linear Schrodinger equation on 
a line. If that soliton is not too big it can be deformed, without prohibitive energy cost, to fit on a unit circle. 
The parameter f3 reflects spontaneous breaking of (ordinary) translation symmetry. Here that breaking is occurring 
through a kind of phase separation. 

Now since non-zero a can be interpreted as magnetic fiux through the ring, we might anticipate, from Faraday's 
law, that as we turn it on, starting from a = 0, our lump of charge will feel a simple torque. (Note that since 
Faraday's law is a formal consequence of the mathematics of gauge potentials, its use does not involve any additional 
hypotheses.) We can also apply "gauge transformations", as in the discussion around Eqn. ([S|). These observations 
are reflected mathematically in the following construction: For any I, we solve 

= l{-^^^-ar^Pl-X\i;l\'^l (21) 

with 

= l{-^^,n-Xm+^^±^i: (22) 

As in the non-interacting ring particle model, the lowest energy is obtained by minimizing Iq + a, for integral Iq. If a 
is not an integer ipi^ (0, t) will be a moving lump, and time translation symmetry will have been spontaneously broken. 
If a is half an odd integer, then its T symmetry is spontaneously broken too. 

This example exhibits several characteristic features of natural r breaking, as discussed in [l| . The lump moves along 
a constant energy trajectory. The parameter /3, which parameterizes an orbit of (ordinary) translation symmetry, 
changes at a constant rate. Although both r and translation symmetry are broken, a combination remains intact. 
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Now let us return to address the conceptual issues alluded to earlier. Our model Hamiltonian was non-local , but we 
required observables to be local. That schizophrenic distinction can be appropriate, since the Hamiltonian might be - 
and, for our rather artificial dynamics, would have to be - carefully engineered, as opposed to being constructed from 
easily implemented, natural observables. Moreover it is not unlikely that the assumption of all-to-all coupling could 
be relaxed, in particular by locating the rings at the nodes of a multidimensional lattice and limiting the couplings 
to a finite range. 

Were we literally considering charged particles confined to a common ring, and treating the electromagnetic field 
dynamically, we would expect radiation from a moving lump of charge. Indeed, the electromagnetic field provides 
modes that couple to all the particles, and in effect provide observers which manifestly violate the framework of Eqn. 
(|12|) . They permit, and will enforce, relaxation to a j/c) state. There are simple variations that can ameliorate this 
issue, e. g. use of multipoles in place of single charges, embedding the system in a cavity to cut off radiation, or simply 
arranging that the motion is slow. A more radical variation, that also would address the unrealistic assumption of 
attraction among the charges, while still obtaining spatial non- uniformity, would be to consider charged particles on 
a ring that form - through repulsion! - a Wigner lattice. 

IV. Imaginary Time Crystals 

In the standard treatment of finite temperature quantum systems using path integral techniques, one considers 
configurations whose arguments involve imaginary values of the time, and imposes imaginary-time periodicity in the 
inverse temperature /? — 1/T. In this set-up the whole action is converted, in effect, into a potential energy: time 
derivatives map onto gradients in imaginary time, which is treated on the same footing as the spatial variables. 

At the level of the action, there is symmetry under translations in imaginary time (iTime). But since iTime 
appears, in this formulation, on the same footing as the spatial variables, it is natural to consider the possibility that 
for appropriate systems the dominant configurations in the path integral are iTime crystals. Suppose that the iTime 
crystal has preferred period A. Then when (3 is an integer multiple of A, the crystal will fit without distortion, but 
otherwise it must be squeezed or stretched, or incorporate defects. Thus we might anticipate periodic behavior in 
thermodynamics quantities in 1/T, with period A, which could serve as an experimental diagnostic. Also, integration 
over the collective coordinate for the broken symmetry makes a contribution to the entropy, even at zero temperature. 
Inspired by the spatial crystal - iTime crystal analogy, one might also be moved to consider the possibility of iTime 
glasses (iGlasses), which would likewise have residual entropy, but no simple order, or iQuasicrystals. 

V. Comments 

1. Spontaneous formation of a time crystal represents the spontaneous emergence of a clock from a time- invariant 
dynamical system. 

2. It is interesting to speculate that a (considerably) more elaborate quantum-mechanical system, whose states 
could be interpreted as collections of qubits, might be engineered to traverse, in its ground configuration, a 
programmed landscape of structured states in Hilbert space over time. 

3. In general, fields or particles in the presence of a time crystal background will be subject to energy-changing 
processes, analogous to momentum-changing Umklapp processes of ordinary crystals. In either case the apparent 
non-conservation is in reality a transfer to the background. In our earlier model, 0{\/N) corrections to the 
background motion arise. 

4. The usual range of questions that arise in connection with any spontaneous ordering, including the nature 
of transitions into or out of the order at finite temperature, critical dimensionality, defects and solitons, and 
low-energy phenomenology, all pose themselves for time crystallization. In addition, there are interesting issues 
around the classification of space-time periodic orderings (roughly speaking, four dimensional crystals). 

5. Semi-macroscopic oscillatory phenomena related in spirit to time crystallization are familiar in the a.c. Josephson 
effect. In that context, however, a voltage difference must be sustained externally. 
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